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ABSTRACT 



The double-lined spectroscopic binary HD 104237 (DX Cha) is part of a complex 
system of some half-dozen nearby young stars. We report a significant change from 
an orbit for the SB2 system derived from 1999-2000 observations. We obtain abun- 
dances from the primary and secondary spectra. The abundance analysis uses both 
detailed spectral synthesis and determinations based on equivalent widths of weak 
absorption lines with Wx typically < 25 mA. Abundances are derived for 25 elements 
in the primary, and 17 elements in the secondary. Apart from lithium and zirconium, 
abundances do not depart significantly from solar. Lithium may be marginally en- 
hanced with respect to the meteoritic value in the primary. It somewhat depleted in 
the secondary. The emission-line spectrum is typical of Herbig Ae stars. We compare 
and contrast the spectra of the HD 104237 primary and two other Herbig Ae stars 
with low V ■ sin(i), HD 101412 and HD 190073. 

Key words: -starsiHerbig Ae -stars:abundances -stars:individual: HD 104237 - 
stars:individual: HD 101412 -stars:individual: HD 190073 



1 INTRODUCTION 

The star HD 104237 (DX Cha) is the brightest of the young 
Herbig Ae stars. Its spatial relationship to several families of 
young stars is nicely illustrated by Feigelson, Lawson, and 
Garmire (2003, see their Figure 1). A detailed description of 
this complex multiple system is available from Grady, et al. 
(2004). One component is very close to the primary so the 
spectrum is of an SB2. 

Bohm, et al. (2004) discuss radial velocities of this sys- 
tem in a "binary approximation," that is, neglecting any 
influence of some half-dozen or more nearby stars. They de- 
rive an orbital period of 19.859 days, and an eccentricity of 
0.665. The systemic velocity was -1-13.94 kms~^. 

A marginal magnetic field detection reported by Do- 
nati, et al. (1997) was not confirmed by Wade, et al. (2007). 
Hubrig, et al. (in preparation) now confirm the presence of 
a weak field. Further discussion of magnetic fields in HD 
104237 and other Herbig Ae stars is given by Hubrig, et al. 
(2007, 2009), and Wade, et al. (2007). 

HD 104237 was among the "24 dusty (pre-)main- 
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sequence stars" whose chemical abundances were deter- 
mined by Acke and Waelkens (2004, henceforth, AW). 
Bohm, et al. (2004) and more recently, Fumel & Bohm (2012, 
henceforth FB) made detailed studies of the pulsation, and 
used the Fe l-Fe ll lines to fix Tcff and \og{g). The present 
abundance study builds on the work of AW and FB. A more 
detailed discussion of the present and earlier work is pre- 
sented in Section 15.11 We do not know of a previous lithium 
abundance determination in a Herbig Ae star. 



2 SPECTRA 

This analysis is based on the spectrum of HD 104237 avail- 
able from the UVESPOP archive (Bagnulo, et al 2003). 
The spectral resolution of the UVES instrument, normally 
80,000, is roughly double that of the material used by AW 
and FB. The signal to noise of the UVESPOP spectra are 
typically between 300 and 500, according to Bagnulo, et al. 
Tests on regions that appeared free of spectral features were 
rarely higher than 200, and were typically between 100 and 
200. 

Additional material from the ESO archives include two 
HARPS spectra, which we designate (a) and (b). HARPS 
spectra have resolution 120000 (Mayor, et al. 2003) but 
the signal-to-noise was significantly lower than that of the 
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Table 1. Observations 



Spectrum 


epoch 


Vh kms-i 


AVh kms-i 


UVESPOP 


23 January 2003 


+36.7 ±0.5 


-20.4 + 1.9 


HARPS(a) 


13 January 2007 


+30.2 ± 1. 


-19.1 + 2.1 


HARPS(b) 


03 May 2010 


-0.6 + 1. 


+37.8 + 1.8 



UVESPOP spectrum. We measured values between 50 and 
60 for HARPS(a). HARPS(b) has S/N 126 in the region of 
the Li I Une, but below 70 in tested regions short of 5200 
A. Therefore, apart from work with the Li I A6707 feature 
fSec. 15.4")) . all synthetic fits and equivalent width measure- 
ments were made on the UVESPOP spectrum. 

Observational epochs are given in Table [1] We also give 
the heliocentric radial velocity of the primary, and the dif- 
ference (Primary minus Secondary) of the radial velocities 
of both components. Errors are based on multiple measure- 
ments of several features. 

The radial velocities of the primary star on 23 January 
2003 and 13 January 2007 are higher than any value mea- 
sured by Bohm, et al., or expected from their orbital param- 
eters: VR{max) ~ +23 kms~^. A significant change in the 
orbit has occurred since the earlier observations of 1999 and 
2000. 
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Figure 1. Synthesized H(5 profile of the composite spectrum 
(gray, or red online) . The observed spectrum (black) shows a par- 
tially filled in, and shifted core. 



The computed spectra for both stars were broadened for 
a gaussian instrumental profile corresponding to a resolving 
power of 80000 for the UVES spectrum and of 120000 for 
the HARPS spectra. 



3 MODEL ATMOSPHERES AND SPECTRAL 
CALCULATIONS 

The abundance analysis was performed by using both spec- 
trum synthesis and equivalent widths. In the first case, 
model atmospheres for HD 104237 (primary P and sec- 
ondary S) were computed with the ATLAS9 code (Castelli 
and Kurucz 2003); in the second case, the T-r relations from 
the ATLAS9 models (Castelli 2010) were used to regenerate 
the depth-dependences based on slightly different opacities 
and partition functions (Cowley, Adelman, and Bord 2003). 
Comparison of the two kinds of calculations showed negligi- 
ble differences in computed features. 



4 ATMOSPHERIC PARAMETERS 

The composite nature of the spectrum of HD 104237 pre- 
cludes the use of broad or intermediate-band photometry to 
fix Tcff and log{g). Therefore we based the parameter deter- 
mination only on the high resolution spectra. 

4.1 Spectral synthesis 

Synthesized spectra of both stars were computed with the 
SYNTHE code (Kurucz 2005). These were added to create 
a synthetic, composite spectrum. This was then overplotted 
with the observed spectrum, and appropriate adjustments 
to the stellar parameters were made to achieve optimum 
agreement. 

Atomic line lists are updated versions of those used by 
Castelli & Hubrig (2004). They are mostly based on the line 
lists available at the Kurucz (2012) website and on the NIST 
(2012) data base. 



4-1.1 The primary 

The derived atmospheric parameters of the primary star 
Tcff ~ 8250 K, and \og{g) = 4.2 were based mostly on a com- 
parison of the observed and computed wings of H5, which is 
a Balmer profile in HD 104237 minimally affected by emis- 
sion (Figure [l|. These parameters were later confirmed by 
the agreement of the Fez and Fell abundances obtained 
from both spectrum synthesis and equivalent widths (see 
below) . 

Because the strong profiles of the primary star were not 
well reproduced after broadening by a rotational profile, we 
added a gaussian macroturbulent velocity. The adopted val- 
ues were v ■ sin(i) of 8 km s~^ for the rotational velocity and 
£,M = 9kms~^ for the gaussian macroturbulent velocity. 

While the weak and medium-strong profiles were well 
predicted by the combined velocity broadenings, the strong 
lines showed an asymmetric profile with a more extended 
violet wing that was impossible to reproduce with our com- 
putations. An example is shown in Figure [2] This kind of 
profile which could indicate an expanding atmosphere, is 
seen in many lines in both the UVES and HARPS spectra. 

The microturbulent velocity value = 2.5kms~^ for 
the primary star was determined both from the equivalent 
widths and the synthetic spectrum methods. 

4-1.2 The secondary 

The spectrum of HD 104237 is contaminated by that of the 
close companion. FB analyzed a spectrum that was obtained 
by subtracting a computed spectrum for the secondary from 
the observed spectrum. 

We have preferred to analyze the observed spectrum as 
it is, and have tried to obtain information also from the sec- 
ondary spectrum. This was not trivial, mostly because no 
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Figure 2. Comparison of the asymmetric observed profile at 
4294.1 A (black line) with the synthetic blend profile due to Till 
4294.094 A and Fel 4294.125 A (gray, red online). 
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Figure 3. Observed (black) and calculated Ca I A6122.21, with 
log(g) = 3.7 (light gray/red online), and log(g) = 4.0 (dark 
gray/blue online). The main feature is due to the primary. Ab- 
sorption from the same line in the secondary star (see text) is to 
the violet. 



usual spectroscopic methods for the parameter determina- 
tions can be applied to the close companion. In fact, neither 
hydrogen lines nor elements with lines in two different ion- 
ization stages are independently observable. We based our 
parameters for the secondary mostly on the Ca I line at 
6122.21 A which is a good luminosity indicator (Cayrel et al. 
1996). The two components for this line are rather well sep- 
arated and free of blends. Unfortunately, the other two lines 
of the Ca I triplet at 6102.73 A and 6161.30 A are blended. 

The composite synthetic spectrum normalized to the 
continuum (PS) was computed according to the relation 
(Lyubimkov & Samedov, 1987) 



F^jPS) _ Fx{P)Rl + F^{S)R's , . 

Fxc{PS) FxciP)Rl + FxciS)Rl' ^' 

We use the subscript 'Ac' to indicate the continuum at a 
wavelength A. The F's for the primary (P) and secondary 
(S) have units of energy per unit area, per sec, per Angstrom. 
The i?'s, are the stellar radii. 

After some experimentation with alternate models and 
different radii ratios, we selected a ratio Rs/Rp = 1 (see 
Section |4.4|) . and the parameters Tcfi =4800K and log(g) 
— 3.7 for the secondary. A rotational velocity vsin(i) = 
12kms~^ was deduced from the comparison of the observed 
and computed profiles. We thank the referee, Barry Smalley, 
who pointed out that the same v sin(i) and macroturbulence 
(8 and 9 kms~^, respectively) would have nearly the same 
combined profile as a vsin(i) of 12 and no macroturbulence. 
This would be the case if the stars were tidally locked. 

A microturbulent velocity ^t = 1.0kms~^ and so- 
lar abundances were assumed. FB used TcS = 4500K, 



log(g) =4.0 and assumed solar abundances and a microtur- 
bulent velocity = 2 km s^^ . 

Figure[3]compares the observed Ca I profile at 6122.21 A 
with the composite profiles computed for Toff — 8250 K, and 
log(g) = 4.2 for the primary and for both T^a = 4800 K, 
and log(p) = 3.7 and T^g = 4500 K, and log(p) = 4.0 for the 
secondary. In the first case the ratio of radii is Rs/Rp ~ 1, 
in the second case it was deduced from the relation 

I§_^Ms_ gp_ 

Rl Mp gs' ^ ' 

for a mass ratio Mp/Ms = 1.29 (Bohm et al. 2004). The 
too-broad wing of the secondary star profile computed for 
log(5r) =4.0 is evident. 

The adopted parameters for the two components (Ta- 
ble. [5} were well suited to reproduce several primary- 
secondary pairs of profiles as Pel at 4950.11, 4969.92, 
5002.793, 5393.17, 5576.09 A, 6065.48, 6252.56, 6335.33 A, 
Cal 6717.68 A, and several others. 

4.2 Equivalent widths 

A closely related approach makes use of measurements in the 
observed composite spectrum (Figure |4)| . Voigt-profile fits 
are used to derive equivalent widths. These values are then 
related to computed equivalent widths from the adopted 
models, assuming solar abundances. 

Equivalent widths were not used to derive atmospheric 
parameters or abundances for the secondary star, with the 
exception of lithium (Sec. 15. 4| ). 

Let xx = L\c{S) /LxciP) be the continuum luminosity 
ratios, where L\c{S) = F\c{S) ■ 4-KRg, and similarly for the 
primary. With the parameters given in Table (2] the fiux 
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Figure 4. Fe I A5391.46, showing Voigt-profile fits to absorp- 
tions from components P and S (weaker, violet component). Lines 
from the secondary can also be seen in the stronger lines, Fe I, 
AA5389.49, and 5393.17 

ratio Fxc{S)/Fxc{P) at A = 5576A is 1/9.9. Because of our 
assumption of equal stellar radii, the ratio is the same for 
the luminosities. 

Directly measured equivalent widths from the compos- 
ite spectrum must be corrected to account for the relative 
contribution of the light of the two stars. Assume that the 
flux in a line of the primary component is diluted only by 
the continuum flux of the secondary component. Then the 
equivalent width of a line of the P-component measured in 
the composite spectrum, Wp^ , is related to the equivalent 
width, Wp in the undiluted spectrum by 

Wp^{l + xx) ■ Wf"" (3) 

Using a similar notation, the (undiluted) equivalent 
width of the secondary component may be obtained from 

Ws = ■ Ws^^. (4) 

Figure 2] illustrates measurements on the composite 
UVESPOP spectrum. It can be seen that the assumption 
that the lines of the two components are separate from one 
another is only approximately true. This assumption is im- 
plicit in Eqs. |3] and U) However, direct comparisons be- 
tween equivalent width and results from complete synthesis 
(including the secondary), are in good agreement, showing 
that the approximation used here are good to a few hun- 
dredths of a dex. At longer wavelengths, the separation of 
the components is greater than that shown in Figure |4l 

Equivalent widths of Fe I and Fe ll lines were measured, 
and corrected as indicated. This is basically an iterative pro- 
cedure. The fraction xx in Eqs. |3] and |4] follows from mea- 
surements of the diluted {Wp^ , and Wg^) and calculations 
of the undiluted {Wp and Ws) equivalent widths. Calcula- 
tions of the latter require models, but initial models are 
available from FB. Adopted models must give similar values 
of Xx from either Eq. [3] or [l] 

If the flux ratios, xx, are known, equivalent widths of 
Fe I and Fe ll lines may be used along with a grid of atmo- 
spheres to determine a relation between T^b and \og{g). 



Table 2. Adopted parameters for primary and secondary atmo- 
spheres 









Cturb 


Primary 


8250±200 


4.2±0.25 


2.5±1 


Secondary 


4800±200 


3.5- < 3.7 


l.Olbl 



We also examined line-by-line abundances for other 
atomic and ionic species with numerous lines, to check that 
the abundances did not drift with excitation potential. This 
provides an independent check on Tcfi 

The overall procedure has the advantage that weak lines 
{Wx ~ 25 mA or less) may be used. These are in principle, 
independent of line broadening mechanisms and the instru- 
mental profile. 

The synthesized spectra were useful for the selection 
of lines that were free from blends, and could be used for 
abundances in the equivalent width method. Indeed, the 
synthetic spectrum automatically includes relevant line ab- 
sorption from the secondary when it overlaps with lines of 
the primary star. Contributions from the secondary are in- 
dicated automatically on the synthesized spectrum, making 
it straightforward to see when a feature is influenced by an 
accidental overlap with an absorption from the secondary 
star. 

The adopted models made use of both spectrum syn- 
thesis and equivalent widths. Their parameters are given in 
Table [21 Uncertainty estimates for the primary are based 
on comparisons of results of calculations based on a grid of 
models with parameters surrounding the adopted values. 

The comparison of the observed profiles of Ca I at 6122 
A with profiles computed for several different choices of Teff 
and logg has lead us to select as best parameters TcS 4800 
K and log{g) <3.7, with an uncertainty of 0.2 dex. While 
there is an uncertainty of about 200K for the temperature, 
the value of the gravity is an upper limit. In fact, values of 
log{g)= 3.6 and 3.5 would have been as well suited to repro- 
duce the observed profile for temperatures differing within 
±200K from 4800K. 

AW used Tcft = 8000K, and log(g) = 4.5. FB's values 
were 8550K, and \og{g) — 3.9, for the primary. 

4.3 Pulsational studies 

Astroseismological measurements are being used to predict 
masses, radii, luminosity, and ages of stars (Mathur, et al. 
2012). Dupret, et al. (2006) attempted to fit the pulsation 
frequencies discussed by Bohm, et al. (2004) to theoretical 
models. They found that the stellar properties that fit the 
observed frequencies were not in good agreement with best 
estimates of these values by other methods. Subsequently, 
Dupret et al. (2007) suggested He accumulation in the par- 
tial ionization zone to account for the observed modes. 

While FB reported extensive pulsational observations 
and analysis, they chose classical spectroscopic methods to 
determine Tea and \og{g) for HD 104237. Those methods 
are basically the same as the ones used in the present work. 
We conclude that astroseismic constraints need to be com- 
plemented by spectroscopic analysis; pulsational properties 
alone are not yet sufficient, at least for HD 104237. 
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4.4 Radii, masses, and consistency 

One approach to the ratio of the radii of the primary and 
secondary of the HD 104237 system is to use the ratio of 
absorption lines, such as those shown in Fig [l] We mea- 
sured equivalent widths of five Fe I hues AA 6065.48, 6252.56, 
6336.82, 6411.65 & 6421.35, to derive the factor xx = 0.155 
of Eqs. |3]and |4l Since that ratio is defined by 



Xx 



LxcjS) 
Lxc{P) 



4 



FxcjS) 
Fxc{P) 



(5) 



we may use the empirically determined xx to obtain the ratio 
of the radii. Note the luminosities and fluxes are wavelength 
dependent. The models of Table [2] provide the surface fluxes 
for each of the stars. That ratio at 6310 A is 0.139, for the 
surface flux of the secondary to that of the primary. With 
these numbers we obtain Rs /Rp — 1.06, or unity within the 
uncertainties discussed in this section and given in Table [S] 
Another way for finding the ratio of radii is the use of 
the Stefan-Boltzmann law. Here, we need the total luminosi- 
ties (integrated over wavelength). Then 



L{S) Rl T^jj[S) 



(6) 



Assuming a luminosity ratio of 10 (Bohm et al. 2004) the 
ratios of radii for T^/ / (P) = 8250 K and T,/ / (S) = 4800 K is 
Rp/Rs = 1.07. 

Yet another approach to finding the ratio of the radii 
is to use the spectroscopically determined surface gravities 
and Eq. [21 Assuming a mass ratio Mp/Ms of 1.29 (Bohm 
et al. 2004), \og{gs) = 3.7, and log(3p) = 4.2, the ratio of 
radii Rp/Rs is 1.6. For Rp/Rs = 1, and the same gravities 
as before, we obtain Mp/Ms =3.15. 

Therefore, while we found a rather good consistency 
between the ratios of radii, luminosities, and effective tem- 
peratures related by the Stefan-Boltzmann law, we obtain a 
manifest inconsistency between the ratios of radii, masses, 
and gravities related by the universal gravitation law. The 
mass ratio obtained by us is in confiict with the value of 1.29 
of Bohm, et al. (2004). However, we must take the uncer- 
tainties into account. 

We note that the spectroscopic gravity for both the pri- 
mary and secondary stars are affected by errors which de- 
pend first of all on the position of the continuum, which 
is rather difficult to draw, in particular above the hydro- 
gen lines. Furthermore, the parameters for the secondary 
are based on too many assumed or estimated quantities to 
allow quantitative analyses of the secondary unaffected by 
some large uncertainties. 

Spectroscopic gravities typically have an uncertainty of 
0.2 (or more) in the log. If we were to assume 4.0 and 3.9 
for the surface gravities of the primary and secondary, then, 
still with unity for the ratio of the radii, we would obtain 
1.26 for the mass ratio. These gravity changes would affect 
the derived abundances by only a few hundredths of a dex. 
Bohm, et al. (2004) assign an uncertainty of 0.02 to their 
mass ratio, so their value could be 1.27. While there is room 
for future adjustments of the fundamental parameters of this 
system, we conclude there is no significant conflict of our 
parameters with the spectroscopic mass ratio. 
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Figure 5. A previously unclassified Fe I line at A4061.09 
(34362.87—58979.82 cm"'^) very closely coincides with a strong 
line of Nd II. Synthesis shows the corresponding feature in HD 
104237 is mostly Fe I. A minor contribution of a few mA, com- 
mensurate with a solar Nd abundance is likely. Labels at the top 
of the plot indicate lines from the primary (gray/red online). The 
decimal fraction of the wavelength is followed by the element (e.g. 
26.01 is Fe ll). The standard SYNTHE plots also contain infor- 
mation on line strengths, and excitation. 



5 ABUNDANCES 

5.1 Previous abundance work 

In this section, we briefly compare and contrast our work 
with that of AW and FB. In the latter study, only the iron 
abundance was determined-independently from Fe I and 
Fe II. 

Both AW and FB study derive abundances from equiv- 
alent widths. Although AW did not account for dilution by 
the secondary spectrum, their abundance results are simi- 
lar to ours. For the well-represented iron spectra (Fe I and 
II ) differences are under 0.1 dex. Differences are less than a 
factor of two, apart from Cu, and two trace species, Nd and 
Sm, where AW report rather large excesses of -1-0.59 and 
-1-0.85 dex. If we average the absolute value of the logarith- 
mic differences, "This study minus AW", the value is 0.13 
dex for the elements through barium. 

FB (see their Table 5) give —4.38 for log(Apc), or 
\og{Fe/H). This becomes -4.42 for \og{Fe/Ntot). This dif- 
fers from our value (Table [4]) by only 0.02 dex. 

Lines from the secondary and weak, previously unclas- 
sified lines could account for the large excesses of Nd and Sm 
reported by AW. An example is shown in Figure[Sl where Fe 
I A4061.09 coincides with Nd II A4061.09. The upper level 
of the Fe I line in question was given by Nave, et al. (1994), 
but is not in Sugar and Corliss (1985). 

FB used abundances from Fe I and Fe ll to fix T^s and 
log(5r). Both FB and AW used excitation to break the TcS- 
log(s') degeneracy. Neither work computed Balmer profiles. 
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Table 3. The effect of microturbulence on abundance 



a(kms-i) 




log(diff) 




4736(74m 


A) 4485(48mA) 


0.0 


0.00 


0.00 


1.0 


0.26 


0.04 


2.0 


0.69 


0.09 


2.5 


0.82 


0.11 


3.0 


0.91 


0.13 


3.5 


0.98 


0.14 



With the help of estimates of the Teff and log(g) of the 
secondary, FB synthesized its spectrum, and subtracted it 
from the observed composite. This method therefore differed 
from the present technique, in which a composite spectrum 
was synthesized and compared directly with the observed 
spectrum. We see no particular advantage or disadvantage 
of subtracting the secondary rather than synthesizing the 
composite spectrum. While FB took the secondary into ac- 
count in their analysis of the primary, they do not report 
abundances for the secondary. 

Given an accurate synthesis of the secondary, FB's 
method assured their equivalent widths were unperturbed 
by blends from the secondary. In the our equivalent width 
method, such blends were avoided with the help of the re- 
sults from the composite synthesis, where contributions from 
the primary and secondary spectrum are indicated on the 
plots (see Figure [5]). 

FB were careful to make use of lines with well-defined 
continua, and accurate oscillator strengths. These con- 
straints, however, reduced the number of usable lines to 8 
Fe I lines and only 2 Fe ll lines. Just 3 of these lines were 
below 30 mA. The remainder were subject to saturation cor- 
rections dependent on the microturbulence, S,t- This was not 
independently determined, but fixed at 2 kms~^, based on 
previous studies. AW used = 3 kms~^. 

Table |3] illustrates the influence of the assumed micro- 
turbulence on abundance for two Fe I lines, A4736 (74 mA) 
and A4485 (48 mA). The entries give the difference in loga- 
rithmic abundance from that obtained with — 0.0. We see 
that for the 74 mA line, the abundance difference using 2.0 
and 2.5 kms"^ is 0.82 - 0.69 = 0.13 dex. For the weaker 
line, the difference is only 0.02 dex. For reference, note that 
in a recent study of Herbig AeBe stars, Folsom, et al. (2012) 
reported values of between 1.3 and 3.7 kms~^. 

We determined independently for each spectrum, obtain- 
ing values ranging up to 3.3 kms~^. Eventually we adopted 
2.5 kms~^, which was more compatible with the spectrum 
synthesis fits. However, our abundances based on equivalent 
widths of weak lines are virtually independent of this choice. 

5.2 Present results 

Equivalent widths, oscillator strengths, and abundances are 
available in the online material for some 375 absorption lines. 
Abundances were obtained from these equivalent widths us- 
ing the Michigan software, usually for for the weak lines 
and are in consistent agreement with best fits from detailed 
spectrum synthesis. 



Table 4. Adopted abundances. Column 2 is a key to oscillator- 
strength references given in the appendix. See text for a descrip- 
tion of the column labeled AW. 



Spec 


gf-Tci 


log(£;//Artot) 


pe 


n 


AW 


Sun 


Li I 


1 


Tabled 










C I 


2 


-3.64 


0.12 


8 


-3.62 


-3.61 


N I 


3 


-4.13 


0.07 


4 


-4.41 


-4.21 


O I 


4 


-3.31 


0.09 


5 


-3.21 


-3.35 


Nal 


5 


-5.64 


0.01 


2 


-5.64 


-5.80 


Mg I 


6 


-4.79 


0.30 


5 


-4.58 


-4.44 


Mg II 


7 


n.u. 






-4.65 




Al I 


8 


-5.66 


0.30 


4 


-5.59 


-5.59 


Si I 


9 


-4.51 


0.44 


7 


-4.47 


-4.52 


Si II 


10 


n.u. 






-4.18 




S I 


11 


-4.70 


0.02 


2 


-4.98 


-4.92 


Ca I 


12 


-5.79 


0.08 


4 


-5.53 


-5.70 


Ca II 


13 


n.u. 






-5.54 




Sc II 


14 


-8.63 


0.20 


5 


-8.87 


-8.98 


Ti I 


15 


-7.05 


0.11 


3 






Ti II 


16 


-7.12 


0.16 


5 


-6.82 




Ti 




-7.08 


0.07 






-7.09 


V I 


17 


-8.06 


0.15 


1 






V II 


18 


-7.77 


0.21 


3 






V 




-7.97 


0.29 




-8.05 


-8.11 


Cr I 


19 


-6.31 


0.25 


7 






Cr II 


20 


-6.26 


0.13 


7 


-6.20 




Cr 




-6.29 


0.05 






-6.40 


Mn I 


21 


-6.63 


0.17 


4 


-6.70 


-6.61 


Mn II 


22 


n.u. 










Fe I 


23 


-4.37 


0.13 


22 


-4.44 




Fe II 


24 


-4.42 


0.17 


17 


-4.42 




Fe 




-4.40 


0.05 






-4.54 


Co I 


25 


-7.26 


0.08 


2 




-7.05 


Ni I 


26 


-5.70 


0.19 


10 


-5.70 


-5.82 


Ni II 


27 


n.u. 




2 






Cu I 


28 


-7.96 


0.25 


3 


-7.57 


-7.84 


Zn I 


29 


-7.58 


0.25 


2 


-7.81 


-7.48 


Sr II 


30 


-9.17: 


0.30; 


2 




-9.17 


Y II 


31 


-9.69 


0.13 


4 


-9.62 


-9.83 


Zr II 


32 


-9.06 


0.06 


4 


-9.09 


-9.46 


Ba II 


33 


-9.61 


0.20 


5 


-9.60 


-9.86 


La II 


34 


-10.72 


0.10 


3 


-10.72 


-10.94 



Abundances are summarized in Table ID Some spectra, 
e.g. Ca 11, were not used for abundances, even though their 
lines were clearly present. This was usually because of blends 
and/or oscillator strengths with low accuracy. Such cases are 
indicated by "n.u." (not used), in the table. Probable errors 
(pe) are standard deviations of abundances from the num- 
ber of lines given in the column labeled 'n'. For n = 2, 
the error is the difference of abundances. For an element, 
e.g. Fe (Fe I and Fe ll), the error is the difference of values 
from the atomic and ionic spectra. The column marked 'AW' 
gives logarithmic abundances computed from AW's paper. 
We used the Anders & Grevesse (1998) solar abundances 
to compute log{El/Ntot) from the bracket values given by 
AW. One exception is Fe, where AW used 7.52 (log(B[) = 12) 
rather than 7.67 given by Anders and Grevesse for the pho- 
tospheric iron abundance. The solar abundances in Column 
7 are from Asplund, et al.(2009). 

Most elements are solar within the errors. We examined 
features near lines of the light lanthanides, Ce-Sm, but they 
were weak and often blended with absorption from the sec- 
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ondary. We see no indication these elements are significantly 
difTerent from solar in abundance. 

A case could be made for slight enhancements of a 
few elements. The case for Zr is the strongest. There are 
numerous lines, modern gf-values (Ljung, et al. 2006) and 
the abundance excess is several times the uncertainty. The 
lithium abundance is discussed below. 



5.3 The secondary 

In addition to Li I, we identified in the secondary lines of 
Nal, Mgl, All, Sil, Cal, Scl, Til, VI, CrI, Mnl, Fel, Col, 
Nil, Cul, Znl, Sri, Ball. A large number of molecular lines, 
especially of C2, CN, CO, CH, MgH contribute to lower the 
continuum. We found several underabundances. Logarithmi- 
cally, relative to solar, they were for Sc [—0.3], Ti [—0.4], V 
[-0.4], Cr [-0.4], Co [-0.7], Ni [-0.3]. Furthermore, in addi- 
tion to Li, we found an overabundance for Ba [+0.5]. All the 
other elements are well predicted by solar abundances. We 
note that these abundances are estimates derived from the 
comparison of the observed UVESPOP spectrum with the 
synthetic spectrum. More work would be needed for an ac- 
curate quantitative analysis of the companion of HD 104237. 



5.4 Lithium 

While lithium is commonly observed in young, cool stars, 
we are unaware of its observation in other Herbig Ae stars. 
Of course, a strong Li I line was noted in HD 104237 by 
FB and others, and correctly attributed to the cooler sec- 
ondary. Bohm, et al. (2004) discuss the lithium line in both 
primary and secondary, but we know of no prior abundance 
determinations from either line in HD 104237. 

Lithium is the only light element that shows a clear 
departure from the solar photospheric abundance. The Li 
I region is shown in Figure (5] The weaker feature marked 
by a vertical line, is Li I in the primary. Note the stronger 
feature, Li I in the secondary, is to the red in HARPS (b). 

An abundance is derived for both primary and sec- 
ondary spectra from equivalent widths as well as synthesis. 
Only abundances from the HARPS (b) spectrum (Figure [6l 
below), where the primary and secondary lines are well sep- 
arated, are adopted. Note that when we apply the correction 
of Eq. |4]to the equivalent width for the secondary, 38.1 mA, 
using xx = 0.171, we get 261 mA, a line no longer weak. 
The abundance, \og{Li/Ntot) = —9.28 was calculated with 
the full hfs structure (cf. Smith, Lambert & Nisson 1998). 
A solar mix of ^Li and '^Li was assumed. 

To estimate errors in the lithium abundances, we use 
HARPS (b) values. Comparing abundance results from from 
synthesis and equivalent widths, we estimate the uncertain- 
ties for the lithium abundances are about 0.1 dex for the 
primary and 0.25 dex for the secondary. 

The lithium abundance for the primary may be slightly 
enhanced with respect to the meteoritic value (Lodders, 
Palme & Gail 2009). There is evidence for some destruction 
of lithium in the cooler secondary. 



I— I 

-d ^ 
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1 1 1 1 1 1 1 1 1 


— 1 1 1 1 1 








UVES ■ 
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6708 
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Wayelength [A] 
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Figure 6. The Li I region, showing synthetic fits. See Section [2] 
for epochs of the UVES, HARPS(a) and HARPS(b) spectra (gray, 
red online) . The vertical line marks the rest wavelength of Li I in 
the primary. Abundances in Table [5] were used for the lower plot. 
Note the different wavelength scales for the upper two, and the 
lower plots. 
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Table 5. Lithium abundances from A6707. Sec text for uncer- 
tainty estimates. 



Primary 




Secondary 




log(Li/Ntot) 


Wx log(Li/Ntot) 


13.6/HARPS(b) 


-8.58 


38.1 -9.28 


synthesis; 






uves 


-8.64 


-8.90 


HARPS(a) 


-8.64 


-8.90 


HARPs(b) 


-8.64 


-9.50 


(solar photosphere:) 


-10.99 




(meteorite:) 


-8.79 





Table 6. Properties of selected emission lines 



FWHM 



Star(HD)/line 


Max / 


w[A] 


[A] 


km s ^ 


Ha 










101412 


3.39 


10.5 


4.42 


202 


104237 


3.55 


13.3 


4.71 


215 


190073 


6.82 


32.2 


5.34 


244 


Na D2 










101412 


1.11 


0.26 


2.36 


120 


104237 


1.45 


1.03 


1.97 


100 


190073 


1.76 


1.38 


1.38 


70 


Fe II 5169 










104237 


1.37 


1.04 


2.17 


126 


190073 


1.77 


1.89 


2.14 


124 


Mg i(bi) 5184 










104237 


1.06 


0.10 


1.48 


86 


190073 


1.45 


1.03 


1.97 


114 


[O] I 6300 










101412 


1.07 


0.125 


1.71 


81 


104237 


1.04 


0.017 


0.75 


36 


190073 


1.11 


0.056 


0.41 


20 



6 THE EMISSION-LINE SPECTRUM 

In this section we compare and contrast the emission-hne 
spectra of HD 104237 with two sharp-lined Herbig Ae stars 
HD 101412 and HD 190073. Spectra and abundances for 
these stars were discussed by Cowley, et al. (2010), and Cow- 
ley & Hubrig (2012). Emission lines increase in prominence 
from HD 101412 to HD 104237 and finally to HD 190073. 
Properties of the more prominent emission lines in these 
three spectra are given in Table (6] The maximum intensity, 
or /-value, and equivalent width, W are with respect to the 
continuum; emission below the continuum line is not taken 
into account. 

Only a few emission lines are seen in HD 101412. These 
include the low Balmer members, the sodium D-lines, and 
He I (Ds or A5876), which is weakly in emission. The strong 
Fe II Multiplet 42, so often in emission in young stars is in 
absorption. Indeed, it is somewhat more strongly in absorp- 
tion than would be expected from the other Fe ll lines, as 
described by Gray and Corbally (2009, see their Nab class). 
Neither the Ca ll H and K lines nor the infrared triplet is 
in emission. 



XI 

<a 



o 




4540 



4550 



4560 
Wavelength [A] 



4570 



4580 



Figure 7. Emission and absorption lines in HD 104237 (center), 
HD 101412 (lower), and HD 190073 (upper). Lines marked with 
filled circles are discussed in the text. 



Both HD 104237 and HD 190073 show the typical emis- 
sion lines of young stars, including the Balmer, Ca ll and 
numerous lines from the first but mostly second spectra of 
iron-group elements. 

Figure [7| illustrates some of the weaker emission lines 
in these stars. The strongest lines in this region are domi- 
nated by Ti II, A4549.60 and 4571.97. Weaker emissions, Fe 
II A4541.56, Cr ll A4558.65, and Ti ll A4563.76 are stronger 
in HD 190073 than in HD 104237. No emissions are evident 
in HD 101412. 



7 DISCUSSION 

Herbig AeBe stars have been investigated as possible progen- 
itors of magnetic Ap/CP stars (Drouine et al. 2004, Folsom, 
et al. 2012). A modest fraction of these stars are magnetic 
(ca. 5-10% Hubrig, et al. 2010, Alecian, et al. 2012), as is 
the case with their main-sequence congeners (Power, et al. 
2008). The chemistry of the Herbig stars is less indicative 
of an evolutionary connection. A few of the Herbigs do in- 
deed have abundance excesses that could connect them to 
classical, magnetic Ap stars (Vick, et al. 2011). However, 
the most common peculiarity that has emerged so far, is a 
resemblance to the A Boo stars. 

In their new detailed study, Folsom, et al. (2012) con- 
cluded that about half of their sample of 20 stars showed A 
Boo characteristics. We had made independent abundance 
studies of two of the Folsom et al. stars, HD 101412 (Cow- 
ley, et al. 2010), and HD 190073 (Cowley & Hubrig 2012). 
The Folsom et al. results for these stars are in excellent 
agreement with ours. These two stars as well as HD 104237 
have relatively sharp spectral lines, an important feature 
for accurate analysis. One star, HD 101412, has mild A Boo 
abundances, while two, HD 190073 and HD 104237 (present 
paper) have solar abundances (apart from Li). Indeed, the 
connection between young stars and the A Boo type has 
recognized for some time (Gray & Corbally 1998). 

The prominence of A Boo types among the Herbig AeBe 
complicates the theoretical understanding of these stars. 
The received mechanism for Ap/CP chemistry, diffusion 
(Michaud 1970), is quite different from that commonly dis- 
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cussed for the A Boo stars, gas-grain separation (cf. discus- 
sion and references in AW). 

HD 104237 has mostly solar abundances, as noted by 
AW. Its excess of lithium is expected for a young star. The 
largest departure from solar reported here is for zirconium, 
whose excess in field stars is usually associated with the 
s-process. If there is an indication of the s-process among 
related elements (Sr, Y, Ba, La), it is marginal. 
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APPENDIX A: OSCILLATOR STRENGTH 
SOURCES 

Keys to oscillator strengths used for the individual spectra 
are given in column 2 of Table |4l A short indication is given 
here to complete citation found in the references. Compila- 
tions of NIST (Kramida et al. 2013), VALD (Kupka, et al. 
1999), and Kurucz (2012) were used for lines not listed in 
modern primary sources. 



APPENDIX B: ONLINE MATERIAL 

Equivalent widths, excitation potentials, oscillator 
strengths, and calculated abundances are available for 
individual lines as online material. 



Table Al. Short References for Table g] 



Key 


Short Ref. 


(1) 


Yan & Drake (1995) 


(2) 


Hibbert, et al. (1993) 


(3) 


Tachiev & Fisher (2002) 


(4) 


NIST/Kramida et al. (2013) 


(5) 


NIST/Kramida et al. (2013) 


(6) 


Kelleher & Podobedova (2008) 


(7) 


NIST/ Fisher, C. F. (2002) 


(8) 


Mendoza et al. (1995) 


(9) 


Nahar & Pradhan (1993) 


(10) 


Kurucz (2012) 


(11) 


VALD Ryabchikova, et al. (1997) 


(12) 


NIST & Aldenius, et al. (2009) 


(13) 


Kurucz (2012) 


(14) 


Lawler & Dakin (1989) 


(15) 


Blackwell, et al. (1982) 


(16) 


Pickering, et al. (2001) 


(17) 


NIST/Ostrovskii & (1958) 


(18) 


NIST & Kurucz (2012) 


(19) 


Sobeck, et al. (2007) 


(20) 


Nilsson, et al. (2006), VALD LS-permitted 


(21) 


VALD/den Hartog et al. (2011) 


(22) 


VALD /LS-permitted 


(23) 


Fuhr&Wiese (2006) 


(24) 


Melendez&Barbuy (2009) 


(25) 


Nitz, et al. (1999) 


(26) 


VALD, Wiclifre&Lawler(1997) 


(27) 


VALD, LS-allowed 


(28) 


NIST/Fuhr & Wiese (2005) 


(29) 


VALD/Kurucz (2012) 


(30) 


Pirronello & Strazzula (1981) 


(31) 


Biemont, et al. (2011) 


(32) 


Ljung, et al. (2006) 


(33) 


NIST/Davidson, et al. (1992) 


(34) 


Lawler, et al. (2001) 



Table Bl. Sample online material for C I. Weak lines used for 
the adopted abundances are marked with an asterisk 

Wavelength EW Log(EW) loggf Chi(eV) Log(C/Ntot) 



4022, 


,844 


9, 


,1 


0, 


,96 


-2, 


.650 


7, 


,480 


-3, 


,71* 


4734, 


,256 


15, 


,8 


1, 


,20 


-2, 


.370 


7, 


.950 


-3, 


,44* 


4770, 


,026 


24, 


,7 


1, 


,39 


-2, 


.440 


7 


.480 


-3, 


,46* 


4932, 


,049 


52, 


,2 


1, 


,72 


-1, 


.660 


7 


.680 


-3, 


,69 


5052, 


, 167 


93, 


,4 


1, 


,97 


-1, 


.300 


7, 


.680 


-3, 


,61 


6010, 


,675 


9, 


,0 


0, 


,96 


-1, 


.940 


8, 


.640 


-3, 
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